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ABSTRACT 

Fatigue-crack  growth  data  for  2024-T3  aluminum  was  used 
to  evaluate  the  accuracy  of  Forman's  equation.   Forman's 
equation  was  found  to  predict  from  one-third  to  twice  the 
observed  cycles  to  a  given  crack  length  when  applied  to 
multiple  cyclic  stress  ratios.   An  alternative  approach 
was  developed,  based  on  the  observed  data,  which  is 
applicable  to  analog  or  digital  fatigue-crack  growth 
prediction.   Data  acquisition  based  on  constant  increment 
of  crack  length  has  been  found  to  yield  inadequate  definition 
of  the  crack-growth  curve  in  the  low-cycle  range  for  the 
successful  application  of  the  method.   An  expression  was 
developed  which  extends  the  application  of  the  crack-closure 
concept  to  the  entire  range  of  cyclic  stress  ratio.   An 
effective  root-mean-square  stress  was  defined  which  shows 
indications  of  good  correlation  with  the  magnitude  of  the 
crack-growth  curves. 
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I.   INTRODUCTION 

Modern  aircraft,  particularly  military  aircraft,  must  be 
designed  for  high  performance  and  long  service  life.   The 
attempt  to  achieve  higher  performance  in  terms  of  increased 
payload,  range,  speed,  or  maneuverability  implies  a  reduc- 
tion in  gross  aircraft  weight,  or  in  the  structural  fraction 
of  that  weight.   This  constraint  results  in  structural 
components  which  carry  higher  loads  per  unit  area  and  places 
greater  demands  on  their  structural  integrity. 

The  design  values  for  the  safe-life  of  non-inspectable 
components  and  the  inspection  intervals  for  other  components 
are  predicated  upon  the  service  life  and  the  predicted  load 
spectrum  specified  by  the  contracting  agency.   The  aircraft 
components  are  sized  to  meet  these  design  values  using 
fatigue  prediction  techniques.   The  estimates  of  remaining 
life  and  maximum  safe  inspection  intervals  of  operational 
aircraft  are  frequently  re-evaluated  in  light  of  the  opera- 
tional load  spectrum  using  similar  methods.   A  factor  of 
safety,  or  factor  of  uncertainty,  is  used  in  these  calcula- 
tions by  which  the  expected  life  of  a  given  structural 
configuration  is  reduced.   This  factor  of  safety  (FS) 
directly  reflects  the  accuracy  of  the  prediction  method  used. 
Large  FS,  indicating  poor  accuracy,  results  in  excessive 
structural  weight  when  compared  to  the  optimum  and  penalizes 
the  aircraft  performance  throughout  the  service  life. 
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Optimistic  estimates  of  the  accuracy  of  a  method  lead 
ultimately  to  costly  retrofits  and  reduced  aircraft 
availability.   The  long  term  economic  efficiency  of  an 
aircraft  is  therefore  dependent  upon  the  accuracy  of  the 
fatigue  prediction  methods  available  to  the  designer  and 
operator;  thus  considerable  motivation  exists  for  the 
improvement  of  fatigue  life  prediction  accuracy. 

The  fail-safe  design  philosophy,  which  admits  the 
existence  of  a  crack  in  the  structure,  is  one  of  the  fati 
prediction  approaches  in  current  use.   Crack  growth  predic- 
tions, in  units  of  cycles  from  an  initial  size  to  its  lengJa 
when  it  precipitates  catastrophic  failure,  are  an  essential 
element  in  the  fatigue  analysis  of  today's  airplanes. 
Inspection  periods  and  rework  intervals  can  be  prescribed 
from  this  information.   This  thesis  reports  an  invest! gat jb& 
of  data  and  theory  used  widely,  if  not  exclusively,  in 
making  crack  growth  predictions  for  aluminum  alloys. 
Specifically,  Hudson's  fatigue  crack  growth  data  [1]  for 
2024-T3  aluminum  alloy  was  examined  in  order  to  establish 
confidence  limits  on  Forman's  equation,  which  appears  to  bji: 
the  most  widely  accepted  crack  growth  prediction  procedure 
used.   Finally,  an  alternate  method  for  expressing  crack 
growth  behavior  is  proposed  to  improve  the  accuracy  of  CESS 
length  predictions. 
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II.   BACKGROUND 

The  fatigue  life  of  aircraft  components  is  difficult 
to  predict  due  to  the  many  factors  influencing  the  structural 
response  to  unsteady  loading.   These  factors  may  be  grouped 
into  four  categories;  material  characteristics,  environmental 
effects,  structural  configuration,  and  service  loading. 

The  fatigue  response  to  unsteady  loads  is  statistical 
in  nature  and  is  dependent  upon  the  time  history  of  the 
loading  and  material  response,  where  material  response  is 
defined  to  be  cycles  to  a  given  crack  length.   The  response 
history  may  be  divided  into  three  stages.   The  first  stage, 
termed  crack  initiation,  is  characterized  by  sub-microscopic 
cracks.   This  stage  has  been  shown  to  be  sensitive  to  a 
variety  of  microscopic  damage  mechanisms.   The  intermediate 
stage,  termed  crack  propagation,  is  characterized  by  the 
relatively  slow  growth  of  a  macroscopic  crack  (a  few  grain 
diameters)  under  the  influence  of  the  ordered  external  load. 
Pinal  failure  occupies  a  small  fraction  of  the  fatigue  life 
of  the  material  and  a  very  high  crack  growth  rate  exists 
during  this  stage. 

A.   MINER'S  RULE  OF  CRACK  INITIATION 

The  obvious  complexity  of  the  problem  has  led  investigators 
to  make  certain  simplifying  assumptions  and/or  approximations. 
Miner  [2]  proposed  a  linear  damage  accumulation  over  the  crack 
initiation  stage  of  the  fatigue  life.   This  proposal  defines 
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damage  as  the  fraction  of  the  initiation  life  In  cycles,  N., 
expended  at  the  I   load  level.   An  increment  of  damage 
would  then  take  place  with  each  cycle  at  the  I   load  level 
according  to  1/N.  or  for  n.  cycles,  the  damage  would  be 
n^/N..   Since  the  accumulation  of  damage  was  assumed  to  be 
linear,  Miner's  Rule  predicts  the  fatigue  life  of  a 
component  as 


L  =  E  n,  (1) 

i   x 


and  damage  as 


ni 

»'i  4  (2> 


At  the  completion  of  the  crack  initiation  stage, 


Z   n,/N.  =  1   .  (3) 

1      1 


Fatigue  life  data  is  presently  available  as  S  vs  N,, 

a     i 

curves  where  S   is  the  alternating  component  of  stress  and 
Nf  is  the  number  of  cycles  to  complete  failure  of  the 
specimen  rather  than  formation  of  a  first  crack.   Application 
of  Miner's  rule  to  random  load  tests  have  yielded  damage 
predictions  at  specimen  failure  considerably  different  from 
unity.   Some  refinements  to  Miner's  approach  to  damage 

prediction  have  taken  the  form  of  weighting  functions  on  each 

mi 
term  or  variable  powers  on  each  term  such  as  D  =  I  (n^/N.) 
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These  modifications  as  well  as  competing  damage  hypotheses 
have  ofttimes  resulted  in  improved  prediction  accuracy  for 
a  given  set  of  data;  however,  the  Miner  summation  of  damage 
continues  to  be  the  most  widely  accepted  due  to  its  simplicity 
and  the  minimal  improvement  in  accuracy  offered  by  its  more 
complicated  alternatives. 

B.   CRACK  PROPAGATION 

After  the  crack  is  formed,  fail-safe  design  practice  of 
crack-growth  predictions  requires  a  knowledge  of  the  crack 
growth  rate  and  the  load  under  which  the  cracked  structure 
would  fail,  both  of  which  are  not  available  from  any  damage 
theory.   Many  investigators,  as  summarized  in  Ref.  3,  have 
attempted  to  correlate  the  crack  growth  per  cycle  to  the 
gross  section  stress  and  the  crack  length,  leading  to 
propagation  models  of  the  form, 

da/dN  =  f(S,a,A)  (4) 

where  2a  is  the  total  crack  length,  S  is  the  gross  section 
stress  and  X  is  a  correction  factor  which  compensates  for 
finite  width,  rivet  forces,  etc.   Paris  and  Erdogan  [4] 
concluded  that  the  stress  intensity  factor  (SIF),  which  is 
a  function  of  the  same  variables,  was  the  parameter  which 
best  correlated  with  the  crack  growth  rate  observed. 
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C.   STRESS  INTENSITY  FACTOR 

Mathematically,  the  stress  Intensity  factors  are  undeter- 
mined coefficients  arising  from  the  solution  of  the  field 
equations  for  any  cracked  geometry.   The  normal  stress  at 
the  crack  tip,  under  moderate  static  loading,  is  given  by 
this  solution  as 


K        ■    K  K 

--,  -  -i  f,;ce)  +  -|  g,,(e)  +4h.,(e)      (5) 


'«     vh  ^r  '     rh  &iru/     vh  "u 


where  K, ,  K~ ,  and  K~  are  the  opening  mode,  in-plane  shear 
mode,  and  out-of-plane  shear  mode  SIP  respectively,  r  is 
the  radial  coordinate  referenced  to  the  crack  tip  and  0  is 
the  angular  coordinate  measured  counter-clockwise  from  the 
advancing  crack  tip.   The  stress  intensity  factors  may 
therefore  be  interpreted  as  parameters  that  reflect  the 
redistribution  of  the  stress  field  in  a  body  resulting  from 
the  introduction  of  a  crack.   The  magnitude  of  the  stress 
intensity  factor  is  dependent  upon  the  geometry  of  the  body 
containing  the  crack,  the  size  and  location  of  the  crack, 
and  distribution  and  magnitude  of  the  external  loads  on  the 
body.   The  configuration  most  readily  treated  is  the  sheet 
specimen  containing  a  central  crack  of  length  2a  under 
Mode  1  loading,  for  which 

K1  =  SAw^a  (6) 
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where  X     is  a  width  correction  factor.   The  criterion  for 
failure  in  the  presence  of  a  crack-like  defect  is  that  crack 
growth  to  failure  (instability)  will  occur  whenever  the 
applied  stresses,  S,  as  given  by  equation  (6)  exceed  some 
critical  condition  specified  by  K,  =  K   .   Paris  and  Ergodan 
[4]  concluded  from  the  preceding  argument  that  this  condition 
may  be  extended  to  fatigue  and  applied  to  the  slow  growth 
phase. 

D.   FORMAN 'S  EQUATION 

Forman  [5]  proposed  a  model  for  crack  propagation  which 
utilized  the  range  of  the  SIF  as  suggested  by  Paris  [6], 


"  =  Kmax  -  Kmln  «) 


to  account  for  the  cyclic  loading  typical  of  the  fatigue 

environment.   This  model  included  the  instability  condition 

for  failure  as  a  singularity  such  that  the  growth  rate 

approaches  infinity  as  K    approaches  K   .   Forman ' s 

J  max  cr 

equation  is  written  as 


da  =    C(AK)n  m 

dN    (l-R)K  r  -  AK  * ° ; 


where  C  and  n  were  assumed  to  be  undetermined  material 
constants  and  R  is  the  ratio  of  minimum  stress  to  maximum 
stress  (loading  ratio). 
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E.   HUDSON'S  DATA 

Hudson  [1]  conducted  a  most  extensive  study  of  the 
material  response  of  two  widely  used  aircraft  alloys, 
2024-T3  and  7075-T6  aluminum,  to  evaluate  Forman's  equation 
when  applied  to  a  large  number  of  maximum  stress  conditions 
and  loading  ratio  combinations.   (The  data  obtained  in 
Hudson's  investigation  was  used  as  the  data  base  for  this 
study. ) 

The  specimens  used  in  Hudson's  study  were  nominally 
12  Inches  in  width,  35  inches  in  length,  and  0.090  inches  in 
thickness.   A  notch  0.10  inches  in  length  by  0.01  inches  in 
width  was  cut  into  the  center  of  each  specimen  by  an  elec- 
trical discharge  machining  technique.   A  reference  grid  of 
0.050  inch  spacing  was  photographically  printed  on  the 
surface  of  each  specimen  to  mark  intervals  in  the  path  of 
the  crack. 

"Fatigue-crack  growth  was  observed  through  10-power 
microscopes  while  illuminating  the  specimen  with  stroboscopic 
light.   The  number  of  cycles  required  to  propagate  the  crack 
to  each  grid  line  was  recorded  so  that  crack-propagation 
rates  could  be  determined.   Approximately  two-thirds  of  the 
crack-propagation  tests  were  conducted  to  failure.   The 
remaining  one-third  were  stopped  before  failure,  and  the 
cracked  specimens  were  used  in  residual-static  strength  tests, 


„1 


National  Aeronautics  and  Space  Administration  Technical 
Note  D-5390,  Effect  of  Stress  Ratio  on  Fatlr.ue-Crack  Growth 
in  7075-T6  and  202'1-T"-;  Aluminum  Alloy  Specimens,  by  C. 
Michael  Hudson,  p.  5,  August  1969. 
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The  results  of  Hudson's  fatigue-crack-growth  tests  on 
the  2024-T3  specimens  are  presented  in  Table  I  (reproduced 
from  Ref .  1) .   Hudson  obtained  the  fatigue-crack  propagation 
rates  graphically  from  the  crack-growth  curves  defined  in 
Table  I.   The  crack-growth  curves  reproduced  from  Ref.  1 
and  shown  as  Figure  1  are  typical.   The  results  were 
presented  as  plots  of  the  common  logarithms  of  the  experi- 
mental value  of  the  crack-propagation  rate  as  functions  of 
the  range  of  the  SIF.   Three  models  were  then  fitted  to  the 
data  by  least-squares  techniques  and  compared  to  the  data. 
Figure  2,  from  Ref.  1,  is  typical  of  the  method  used  for  the 
presentation  of  data  and  the  evaluation  of  the  quality  of 
fit  of  empirical  equations.   Hudson  concluded  that  Forman's 
equation  produced  an  excellent  fit  for  both  materials  whereas 
the  other  equations  evaluated  "...produced  good  correlation 
with  the  test  data  except  at  the  high  growth  rates  for  the 
7075-T6  alloy"  on  the  basis  of  graphical  representations  of 
which  Figure  2  is  typical.   This  conclusion  at  first  appears 
to  be  justified  on  the  basis  of  plots  such  as  Figure  2, 
however,  the  fatigue  analyst  is  concerned  with  the  accuracy 
with  which  a  crack-propagation  model  predicts  the  number  of 
cycles  required  for  the  crack  to  attain  a  given  length.   This 
is  the  basis  upon  which  the  fit  should  be  evaluated. 
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III.   DEVELOPMENT 

A.   EVALUATION  OF  FORMAN'S  EQUATION 

Hudson  used  the  tangent  form  of  the  width  correction 
factor, 


Xw  = 


W  .  iTa 
— —  tan  -rr 
■rra      W 


(9) 


in  evaluating  the  material  constants,  C  and  n,  in  Forman's 
equation.   The  differential  equation  resulting  from  the  use 
of  this  form  of  width  correction  in  Forman's  equation  may 
be  integrated  (Appendix  A)  using  the  small  angle  approximation, 
tan  6  ^  9  .   This  approximation  is  accurate  to  within  three 
percent  for  crack  lengths  of  less  than  0.8  inches  for 
Hudson's  specimens. 

The  predicted  interval  of  cycles  between  crack  lengths, 
AN  .  was  calculated  using  Hudson's  fit  of  Forman's  equation 
and  the  integrated  form  of  Forman's  equation  becomes 


AN,  =  2C\ 
P     1 


where 


"   2=H  K     S   \/a 

2  (    cr    max  V  o^ 

ao    Cn~=2     H^3    ] 


(1-R) 


2-n 


K 


2  ,   cr 

a     (?T=2 


S   Va 

max  v   v 

n-3   . 
(10) 


C[(1-R)S    ] 
max 


n 
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and  AN  is  the  predicted  number  of  cycles  required  to 
propagate  at  half-crack  of  length  a  to  length  a.   The  ratio 
of  the  predicted  to  the  observed  cycle  increment  is  presented 
in  Table  II,  where  gross  discrepancies  exist  for  almost  every 
case. 

The  load  cases  were  also  fitted  individually  in  order  to 
determine  the  cause  of  the  poor  predictions  reported  in 
Table  II.   The  expression  for  AK  used  in  this  analysis  was 


AK  =  (S    -S  .  )  A  -/-^  (11) 

max   min   wV  Z 


where  Brown  and  Srawley's  [7]  finite  width  correction  for  a 
centrally  cracked  sheet, 


Xw  =  1.77  +  -277(^)  -  .5K^)2  +  2.7(^)3,     (12) 


was  chosen  on  the  basis  of  its  reported  accuracy.   The  test 
data  for  each  load  case  were  represented  by  a  power  series 
of  the  form 


a  =  p-jN  1  +   p2N  2  +  ...  +  PmN  m  (13) 


which  yields  a  very  good  approximation  of  the  observed  data 
as  indicated  by  Figure  3.   An  analytic  crack-growth  rate 
expression  was  obtained  by  performing  a  differentiation  of 
equation  (13)  with  respect  to  cycles  leading  to 
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§r  =  Pi6/1""  +  p?^^2"1  +  •••  +  Vr/111"1     (12|) 


Equations  (11)  and  (13)  were  applied  to  each  load  case  and 
a  least-squares  fit  of  Forman's  equation  to  the  data  was 
performed.   Figure  4  is  typical  of  the  individual  fits  thus 
obtained. 

A  regular,  periodic  variation  of  the  calculated  crack- 
growth  rate  about  the  fitted  line  was  observed  which  is 
believed  to  indicate  that  Forman's  equation  may  be  an 
incorrect  or  incomplete  functional  form.   These  periodic 
variations  appear  as  scatter  when  Forman's  equation  is  called 
to  fit  multiple  load  cases.   An  apparent  lower  bound  of 
crack-growth  rate  may  also  be  observed  in  Figure  4.   The 
failure  of  Forman's  equation  to  represent  this  threshold 
level  probably  accounts  for  the  largest  portion  of  prediction 
error.   Just  as  for  the  integration  performed  for  Hudson's 
fit,  consistent  prediction  errors  were  noted  for  the  fits 
of  the  individual  load  cases  but  of  smaller  magnitude.   The 
periodic  variation  of  the  data  about  the  fitted  line  and  the 
apparent  lower  threshold  of  crack-growth  rate  led  to  the 
development  of  a  new  material  response  model.   Equation  (13) 
was  not  considered  to  be  a  practical  alternative  to  Forman's 
equation  in  view  of  the  explicit  appearance  of  the  variable 
N,  which  is  unknown  in  the  practical  application  of  any 
prediction  model.   Series  reversion  was  likewise  discounted 
due  to  non-integer  exponents  resulting  from  the  fit. 
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B.   MATERIAL  RESPONSE  FUNCTION 

A  material  response  function  was  developed  in  which  the 
response  was  assumed  to  be  a  piecewise  continuous  function 
of  time,  loading  conditions,  and  crack  geometry.   The 
function  was  also  assumed  to  possess  piecewise  continuous 
time  derivatives.   This  assumption  was  based  on  the  satis- 
factory representation  of  the  data  achieved  by  the  power 
series  approximation  of  the  fatigue  data,  equation  (13) s 
which  has  an  infinite  number  of  continuous  time  derivatives. 

Previous  investigators  have  assumed  that  the  fatigue- 
crack  propagation  rate,  da/dN,  has  the  functional  form 

da/dN  =  f(S(N),a(N),X(N))  (15) 

where  the  variables  have  the  meanings  previously  defined  and 
are  functions  of  time  or,  alternatively,  cycles.   Although 
the  relationship  between  the  independent  variables  N  and  t 
is  linear,  as  expressed  by 

N  =  ut  (16) 

where  to  is  the  circular  frequency  of  loading,  a  distinction 
in  the  growth  rates  da/dt  and  da/dN  exists.   The  crack-growth 
rate  da/dt  is  an  instantaneous  rate  which  must  be  integrated 
over  the  period  of  a  cycle  in  order  to  obtain  the  time 
average  rate,  da/dN. 
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Paris  [6]  concluded  that  the  crack  growth  rate  was  best 
represented  as  a  function  of  the  range  SIF  as  given  by 
equation  (11).   The  desired  function  must  predict  crack 
growth  rate  under  constant  amplitude  and  random  loading  with 
equal  precision  since  constant  amplitude  loading  is  a  special 
case  of  random  loading.   The  stress  range  is  therefore  a 
constant  for  either  load  case  taken  on  a  cycle-by-cycle  basis. 
This  argument  applies  equally  well  to  any  other  choice  of 
the  characteristic  stress  to  be  used  in  calculating  the  SIF. 
The  expression  for  crack-growth  rate  could  then  be  represented 
equally  well  by  the  functional  form 


§§=  f1(a(N)3Xw(N))  (17) 


or 


f§  =  f2(x(N))  (18) 


where 


x  =  Aw(a,W)-\/a(N) 


since  X  is  a  function  of  crack  geometry  only. 

Introduction  of  the  requirement  for  higher  order  deri- 
vatives can  be  satisfied  by  a  linear  summation  of  the  assumed 
form 


cdx  +   cd_j£+  +cd_21  (19) 

Cl    dN        c2     ,..2         •••Ck,k  KXJ} 

dN  dx 


2 'I 


where  the  c1  =  ci(S).   Equation  (19)  is  a  linear  differential 
equation  In  N  with  constant  coefficients  and  thus  has  a 
solution  of  exponential  functions.   Paris  [6]  noted  that  the 
crack  growth  rate  is  a  monotonically  increasing  function 
of  the  stress  intensity  factor.   The  crack  length,  a, 
increases  monotonically  with  time  and  the  width  correction, 
X,  increases  monotonically  with  crack  length,  therefore  the 
rate  of  change  of  the  variable',  x,  with  respect  to  cycles 
(time)  must  increase  monotonically  with  x.   This  type  of 
response  is  typical  of  positive  feedback  systems,  which 
implies  that  the  roots  of  the  characteristic  equation 
obtained  from  equation  (19)  are  real.   A  mathematical  model 
of  the  desired  material  response  function  is  therefore  the 
Prony  series 


x  =  Z      b2i-l  exp[b2iN]  (20) 


which  satisfies  all  requirements  heretofore  set  forth. 

This  form  of  the  material  response  function  may  be  recast 
to  the  differential 


,k  ,  ,  ,2  ,k-l 

dx_lr  dx    ,  d  x    ,  ,  d  n 


dNk        ck    L    -1    dN         -2    dN2         ' • •         -k-1   Nk-1 


(21) 


when  the  functional  form  of  the  c.  is  determined.   This 
differential  equation  would  find  utility  if  programmed  for 
analog  or  hybrid  computation  using  non-linear  elements 
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eliminating  the  excessive  digital  computer  time  requirements 
typical  of  random  loading  fatigue  calculations;  however,  it 
is  beyond  the  scope  of  this  thesis.   This  approach  to  fatigue 
prediction  calculations  would  also  allow  the  use  of  discretized 
flight  load  data  without  intermediate  processing. 

The  data  sampling  frequency  used  by  Hudson,  based  on 
equal  increments  of  crack  length,  provided  insufficient  data 
in  the  low  cycle  range  of  the  specimen  lives.   The  crack 
growth  rate  progresses  from  a  nearly  constant  value  at  small 
crack  lengths  to  much  larger  values  within  a  small  increment 
of  crack  length  requiring  a  more  frequent  data  sampling 
frequency  in  order  to  achieve  a  successful  analysis  of  the 
exponential  form  of  the  fatigue  response  function  through 
direct  analysis. 

C.   EFFECTIVE  STRESS 

The  extension  of  equation  (20)  to  account  for  random  load 
sequences  requires  that  the  constants  be  a  function  of  the 
loading  conditions.   The  representation  of  the  piecewise 
continuous  response  function,  x'(t),  as  the  cycle  average 
function  x(N)  leads  to  the  representation  of  the  piecewise 
continuous  crack-growth  rate,  da(t)/dt,  as  the  time  average 
crack-growth  rate,  da(N)/dN.   Consistency  requires  that  the 
loading  also  be  represented  as  a  time  average  variable. 
Paris  [16]  suggested  the  utility  of  the  Root-Mean-Square 
(rms)  gross  section  stress,  as  defined  as 
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S™s2  "  k     {     s2(t)  «  <22> 


where    S(t)  =  Sm  +  Sa  sin  ut 

and  Sm  is  the  mean  gross  section  stress  during  the  cycle, 
Sa  is  the  amplitude  of  the  gross  section  varying  stress 
during  the  cycle,  w  is  the  circular  frequency  and  T  is  the 
period  of  the  cycle. 

The  entire  magnitude  of  the  stress  described  above  does 
not  contribute  to  crack  propagation.   Elber  [8]  has  noted 
that  plastic  yielding  at  the  crack  tip  results  in  a  compres- 
sive residual  stress  upon  unloading.   This  compressive  stress 
at  the  crack  tip  delays  the  opening  of  the  crack.   Elber 
defined  the  crack  closure  parameter,  U,  as 


U  =  3max  "  g0P  (23) 

max    min 


where  S   is  the  gross  section  stress  at  which  the  crack 
opens.   He  found  that  U  seemed  to  be  a  linear  function  of 
loading  ratio  for  2024-T3  aluminum  as  represented  by 

U  =  0.4R  +  0.5  (2^) 

Since  no  crack  propagation  occurs  while  the  crack  is 
closed,  an  effective  rms  gross  section  stress  may  be  defined 
as 
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Srms2  "  I     Se2  f  ^ 


where       S  =  Sm  +  S  sin  ?J*  _  s^ 
e    m    a     T     op 


subject  to  Sg  >  0  . 

The  nature  of  the  sinusoidal  function  is  such  that  the  limits 
of  integration  can  be  shifted  to  write  (25)  more  conveniently 
in  terms  of  the  damage  effective  portion  of  the  cycle  as 

illustrated  in  Figure  5.   This  leads  to 

1  +_f_  1  _f_ 

f  2T  ¥  2T 

s    2  =  f      s2  &-      +  f  o4r   (26) 

rms     J  e   T     J  T 

'  1  _f_  1  _f_ 

¥  2t  V     2t 


where  f/T  is  the  fraction  of  the  cycle  during  which  the 
crack  is  open  and  is  defined  by  (See  Appendix  B) 


f/  g  cos  X(1-2U)  .      . 


Elber's  equation  for  U(R)  leads  to  effective  cycle 

fractions  which  are  undefined  for  load  ratios  less  than  -1.25 

This  unnatural  behavior  in  the  neighborhood  of  a  subset  of 

Hudson's  data  motivated  a  search  for  a  more  satisfactory 

representation  of  the  crack  closure  phenomenon.   Figure  7 

is  a  schematic  of  the  loading  conditions  for  increasingly 

negative  load  ratios.   Reasoning  from  this  figure  indicates 

that,  in  the  limit,  as  S    approaches  zero  (R  approaches 

max 
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negative  infinity) ,  the  fraction  of  the  period  over  which 
the  crack  is  open  approaches  zero.   This  physical  constraint 
is  satisfied  by  an  equation  of  the  form 


T  ~  B^R  (28) 


where  A  and  B  were  assumed  to  be  undetermined  material 
constants.   Substitution  of  equation  (28)  into  the  defini- 
tion of  f/T,  equation  (27),  leads  to 


U  =  |  [1  -  cos  tt  (^Ar)]  (29) 


Elber's  data  [8]  was  fitted  to  the  function  for  the 
effective  cycle  fraction,  equation  (28),  as  shown  in  Figure 
6.   The  form  satisfies  the  physical  constraint . and  yields 
a  smaller  x   error  than  equation  (24)  when  compared  in 
Figure  7.   Figure  7  is  a  comparison  of  the  functional 
representation  of  U(R)  as  developed  in  this  study  with  that 
of  Elber.   Both  forms  appear  to  be  incomplete  in  that  neither 
form  considers  the  effects  of  absolute  load  level  or  crack 
geometry  on  the  crack  closure  phenomenon.   S     was 

I  III t>  q 

calculated  on  the  basis  of  equations  (25)  and  (28). 


D.   HYPERBOLIC  REPRESENTATION  OF  FATIGUE -CRACK  GROWTH  DATA 
On  the  basis  of  the  arguments  presented  with  regard  to 
the  exponential  form  of  the  material  response  function  and 
the  observed  similarity  of  the  exponential  and  hyperbolic 


29 


curve  forms,  Hudson's  fatigue-crack  growth  data  was  fitted 
to  the  general  form  of  the  second  order  equation. 


Ax2  +  BxN  +  CN2  +  Dx  +  EN  =   1  (30) 


Solving  for  N,  Table  III  is  a  listing  of  the  coefficients  of 
the  equation 


N  =  d-jX  +  d2  ±  , y&Tx2  +  d^x  +  d5    ,  (3D 

ordered  with  decreasing  S    ,  which  resulted  from  the  fit 

rms 

of  equation  (30) . 

Figure  8  indicates  the  quality  typical  of  these  coeffi- 
cients.  The  variation  in  sign  which  occurs  in  Table  III 
reflects  both  the  translation  of  the  x  axis  (or  the  N 
ordinate  of  the  center  of  the  fitted  hyperbola)  and  the 
frequent  occurence  of  parabolic,  rather  than  hyperbolic, 
coefficients.   The  parabolic  coefficients  (sign  of  A  =  sign 
of  C  in  equation  (30))  resulted  from  the  paucity  of  data  in 
the  low-cycle  range  of  the  specimen  fatigue  life.   The 
appearance  of  parabolic  coefficients  is  unacceptable  from 
a  physical  standpoint  as  the  parabolic  form  implies  that 
cycles  to  a  given  crack  length  is  not  a  single  valued 
function,  as  shown  schematically  in  Figure  9a.   The  trans- 
lation of  the  center  of  the  assumed  hyperbola  was  an 
unavoidable  result  of  the  presence  of  the  original  machined 
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notch.   Figure  9b  Is  a  schematic  representation  of  the 
effect  of  load  level  and  the  presence  of  the  artificial 
notch  on  the  translation  of  the  x  axis. 

E.   NORMALIZATION  OF  THE  DATA 

A  unified  fatigue  prediction  method  may  be  formulated 

by  combining  the  capabilities  of  the  damage  accumulation 

method  in  estimating  crack  initiation  with  the  fracture 

mechanics  approach  applied  over  the  crack  propagation  and 

failure  stages.   This  usage  demands  the  normalization  of 

all  load  cases  to  a  common  reference  crack  length  which 

should  be  the  minimum  detectable  crack  length  defining 

initiation.   This  normalization  was  accomplished  for  each 

(hyperbolic)  load  case  by  estimating  the  cyclic  interval 

between  Hudson's  reference  point,  typically  0.10  inch,  and 

the  assumed  initiation  points,  0.003  inch.   The  estimate  was 

made  on  the  basis  of  equation  (31)  using  the  d.  previously 

calculated  and  was  then  applied  as  an  increment  to  the 

observed  cycles  to  a  given  crack  length  for  an  unnotched 

specimen.   The  modified  data  sets  were  again  fit  to  equation 

(30)  yielding  estimated  coefficients,  d  .*  ,  of  the  unnotched 

sheet  which  are  hypothesized  to  correlate  with  S 

rmse 

Figures  10  through  Ik    indicate  the  degree  of  this  correlation 
and  suggest  that  the  appropriate  functional  relation  may  be 

V  "  »i  W"1  (32) 

e 
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An  analytic  expression  for  the  crack-growth  rate,  da/dN, 
may  be  obtained  by  the  chain  differentiation 


da  _  da  dx  _   1    1  fr>o\ 

dN    dx  dN  "  ,dxv  rdNx  UJ; 


where   g.-±-(*  +  af|)  (3*0 


|  (2d  x  +  d^) 
and     jg-  =  dl  -  r  (35) 

rd_x  +  d^x  +  dj- 


11  "7 


This  expression  for  da/dN  approaches  the  constant  (threshold) 
crack-growth  rate  observed  in  Figure  *i  at  small  crack  lengths 
and  satisfies  the  instability  condition  given  that  equation 
(3*0  has  a  zero. 

The  zero  of  equation  (3*0  yields  x.    and  leads  to 


K.  „  =  S    x.  (36) 

ins    max  ins 


where  the  instability  condition  has  been  assumed  to  take 

place  at  the  instant  in  the  last  cycle  at  which  the  load 

has  attained  its  maximum  value.   Figure  15  depicts  the 

variation  of  the  ratio  K.   /K   .   This  apparent  reduction  of 

ins   c  3? 

the  level  of  the  stress  intensity  factor  for  unstable  crack 
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propagation  is  attributed  to  dynamic  effects  on  propagation 
and  has  the  effect  of  minimizing  the  influence  of  the 
singularity  condition  in  Forman's  equation. 

The  data  fits  to  equation  (30)  yielded  no  real  solution 
to  equation  (3*0  for  some  load  cases.   These  load  cases 
appeared  to  be  the  dynamic  tests  which  were  terminated  well 
before  failure  and  therefore  contained  Insufficient 
Information  regarding  the  shape  of  the  response  function  in 
the  high-cycle  region  to  give  the  correct  physical  response. 
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IV.   CONCLUSIONS 

Hudson's  [1]  fatigue-crack  growth  data  for  2024-T3 
aluminum  was  used  to  evaluate  the  accuracy  of  Forman's 
equation.   An  alternative  approach  to  the  representation 
of  the  propagation  process  was  developed  and  analyzed.   An 
expression  was  developed  which  extends  Elber's  crack-closure 
concept  to  the  possible  load  ratio  range.   The  following 
conclusions  can  be  drawn  from  this  study. 

1.  Forman's  equation  does  not  appear  to  satisfactorily 
predict  crack  propagation  in  2024-T3  when  evaluated  on  the 
basis  of  the  observed  data. 

2.  Elber's  crack-closure  concept  was  successfully 
extended  to  the  entire  load  ratio  range  in  2024-T3  by  the 
relation 

U  -  Sln   (2.l69-R} 

3.  The  development  of  a  crack-growth  rate  expression 
based  on  empirical  analysis  of  the  observed  data  met  with 
limited  success  due  to  the  coarse  data  resolution  in  the 
low-cycle  range. 

4.  The  stress  intensity  factor  at  crack  instability 
under  dynamic  loading  conditions  was  observed  to  decrease 
with  decreasing  values  of  effective  rms  gross-section 
stress,  asymptotically  approaching  a  value  of  approximately 
one-third  the  value  measured  under  static  loading. 
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IV.   RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

The  material  response  of  2  02  4-T3  aluminum  has  been  shown 
to  be  capable  of  representation  by  a  continuous  response 
function.   The  accumulation  of  a  number  of  test  sets  with 
more  frequent  data  sampling  in  the  low-cycle  phase  of  the 
crack  propagation  life  would  allow  a  determination  of  a 
material  response  model  for  which  statistical  reliability 
limits  could  be  determined.   Any  such  model  must  be  based 
upon  the  observed  data,  as  compared  to  derived  data,  for 
the  reliability  limits  to  be  meaningful. 

The  development  of  the  exponential  response  model  will 
allow  the  determination  of  the  governing  differential  equation 
for  crack  propagation  under  dynamic  loading  by  the  application 
of  parameter  identification  techniques.   Such  a  development 
would  allow  the  calculation  of  the  expected  crack  length 
under  random  loading  by  the  use  of  analog,  rather  than 
digital,  computation  devices. 
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APPENDIX  A.   INTEGRATION  OF  FORMAN'S  EQUATION 


Forman's  equation  is  written 


da  _   C(AK)n  ,,   . 

dN  ~  (l-R)K   -  AK  KR   ±} 


where  the   stress   intensity   factor,   K,    is   defined  as 


K  =■   S\-\fta  (A-2) 


and  the  range  of  the  stress  intensity  factor  is  written 


AK  =  ASX^N/Jra"  (A-3) 


The  tangent  form  of  the  width  correction  factor, 


*w  -\/?  tan  if  (A~^ 

w    V  -rra      W 


may  be  used,  allowing  the  small  angle  approximation.   This 
approximation,  tan  0=9,  yields  errors  of  less  than  three 
percent  for  ratios  of  crack  half-length  to  width  of  less 
than  1/15. 

Limiting  the  integration  to  this  range  allows  AK  to  be 
expressed  as 

AK  =  AS\/a~  (A-5) 
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The  cyclic  stress  ratio,  R,  is  defined 


R  = 


max 
5min 


(A-6) 


Therefore   AS  =  (l-R)S 


max 


(A-7) 


Forman's  equation  can  be  rewritten  as 


da   ^^  Smaxn  ^ 

dN    (iI^c7rwI7r) 


(A-8) 


for  the  limitations  imposed. 

Treating  this  expression  as  a  differential  as  proposed 
by  Porman  [5] 


(1-RMK   -  S    a1/2) 
dN  =  cr  -  P»; da 

C(l-R)nS   na1/2 
max 


(A-9) 


which  may  be  integrated  between  crack  half-lengths  a  and  a 


AN   = 


(1-R) 


P    C(l-R)nS 


-  V- 

n    J 


I 

max     o 


K   -  S  n  a 
cr   max 


1/2 


a 


HTT 


da 


(A-10) 


Carrying  out  the  indicated  integration  yields 


ANp    ■ 

(1-R) 

r              2~n 

K        a   2 
cr 

-      2-n 

3-n     _ 

S          a2 
max 

3-n            J 

a 

C(l-R)nS 

n 

1  V 

a 

(A-ll) 
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Simplification  of  equation  (A-ll)  leads  to 


r   2-n 


AN   =  2C, 
P      1 


n-2 


max 


^Fo 


n-3 


2-n 


-  a 


n-2 


max 


A^" 


n-3 
(A-12) 


where 


Cl  = 


(1-R) 


C(l-R)nS   n 

max 


(A-13) 


and  AN  is  the  cyclic  increment  predicted  to  elapse  during 

the  extension  of  a  crack  of  half-length  a  to  a  half-length 

of  a  when  subjected  to  constant  amplitude  loading  conditions 

R  and  Smax- 
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APPENDIX  B.   THE  DERIVATION  OP  THE  RELATIONSHIP  BETWEEN 
THE  DAMAGE  EFFECTIVE  CYCLE  FRACTION  AND 
THE  CRACK-CLOSURE  PARAMETER 


Crack  propagation  occurs  only  when  the  fatigue  crack 
is  open,  allowing  the  definition  of  a  damage  effective 
portion  of  a  cycle  (period) .   This  fraction  of  the  period 
is  bounded  by  the  limits  t./T  and  tp/T  representing  the 
instants  of  crack  opening  and  crack  closing,  respectively. 
The  damage  effective  fraction  is  then 


f/T  =  t2/T  -  t^/T  (B-l) 


These  limits  arc-  symmetric  about  the  maximum  load  instant 
for  constant  amplitude  loading.   Assuming  a  loading  function 
of  the  form 


S(t/T)  =  Sm  +  So  sin  (27Tt/T)  (B-2) 

m    a 


the  limits  of  the  damage  effective  fraction  of  the  cycle 
may  be  written 


t-j/T  =  1/4  -  f/(2T)  (B-3a) 


and 


t2/T  =  1/H  +  f/(2T)  (B-3b) 
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The  crack  opens  during  the  loading  portion  of  the  period 
when  the  effective  gross  section  stress  is  zero.   The  crack- 
closure  concept  modifies  this  stress  by  the  inclusion  of 
the  residual  stresses  resulting  from  local  tensile  yielding 
at  the  crack  tip.   Defining  S   to  be  the  gross  section 
stress  at  which  the  crack  opens  leads  to 


So  (t/T)  -  S  .-  S_  +  SQ  sin  (2-irt/T)         (B-4) 
e  m    op    a 


The  instant  of  crack  opening  is  therefore  defined  by 


So  (t,/T)  =  Sm   -  Snn   +   SQ  sin  (2irt,/T)  =    0    (B-5) 
c    j.        m     op     a.  j. 


or 


S    —  S 

(^S ~]    =    Sln  l^2^^/4  -  f/2T)]  (B-6) 

a 


Trigonometric  reduction  yields 


S    -  S 
(_o2^ m}  =  cos  (7rf/T)  (B_7) 


and  the  damage  effective  fraction  of  the  cycle  is  defined 
In  terms  of  the  loading  and  closure  parameters 


S    -  S 

f/T  A  cos"1/!!  (-2E. 01)  (B-8) 

a 


'10 


Rewriting  Equation  (B-8)  in  terms  of  S  „   and  S  . 

max      mln 


yields 


f/T  = 


cos  /it 


S    +  S 

e,     /  max min n 

op  "  [  2      J 

max mm; 


(B-9) 


which  simplifies  to 


f/T  =  cos   /it 


1  - 


2(S     -  S   ) 

max    op 

S    —  S 
max    min 


(B-10) 


This  becomes 


f/T  =  cos~1/tt  [1  -  2U] 


(B-ll) 


with  Elber's  [8]  definition  of  the  crack-closure  parameter 
U, 


s    —  s 

TT  .   max op 

S    —  S 
max    min 


(B-12) 


Solving  Equation  (B-12)  for  U  yields 


U  =  i  (1  -  COS  TTf/T) 


(B-13) 


Therefore,  the  relation  between  f/T  and  U  is 


U  =  sin   (^r) 


(B-1'0 


41 


and 


f/T  =   2/tt   sin"1    (U)2 


(B-15) 
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TABLE  II.   Analysis  of  the  prediction  accuracy  of  Forman's 
equation  for  small  cracks  in  2024-T3  aluminum. 
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Figure  2.   Hudson's  comparison  of  fatigue-crack  growth  rate 
with  Forman's  equation  (From  Ref.  1). 
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effective  W.S   gross-section  stress 
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